FLUCTUATIONS OF EIGENVALUES OF RANDOM NORMAL MATRICES 

YACIN AMEUR, HAKAN HEDENMALM, AND NIKOLAI MAKAROV 

Abstract. In this note, we consider a fairly general potential in the plane and the corresponding 

Boltzmann-Gibbs distribution of eigenvalues of random normal matrices. As the order of the 

matrices tends to infinity, the eigenvalues condensate on a certain compact subset of the plane - the 

r*0 | "droplet". We give two proofs for the Gaussian field convergence of fluctuations of linear statistics 

of eigenvalues of random normal matrices in the interior of the droplet. We also discuss various 

^2 ' ramifications of this result. 
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1. Notation, preliminaries and the main result 

Random normal matrix ensembles. Let a weight function (or potential) Q : C — > Rbe fixed. We 
Oh ■ assume throughout that Q is C°° on C (sometimes excepting a finite set where the value may be 

C"| , +°°) and that there are positive numbers C and p such that 

(1.1) Q(z)>plog|z| 2 , |z|>C. 

Let 9J„ be the space of all normal n x n matrices M (i.e., such that M*M — MM*) with metric 
induced from the standard metric on the space C" of all n x n matrices. Write M = UDU* where 
U unitary, i.e. of class lt„, and D = diag(A,) e C". 

It is well-known 11151 , 1118 1 that the Riemannian volume form on 5R„ is given by AM n := 
dil„ I Vn(Ai, • • • , A„)| 2 d 2 Ai • • • d 2 A„, where dil n is the normalized U„ -invariant measure on M„/T, 
and V n is the Vandermonde determinant 



^1" 



r^ 



V t ,(A 1 ,...,A„) = Yl(A j -A k ). 

o ; i<k 

00 . 

We introduce another parameter m>\ and consider the probability measure (on s Jt„) 






dP m ,„(M) = -!- e -'" tra «Q(M) dMn/ 



w^here C„, „ is the normalizing constant making the total mass equal to one. 

In random matrix theory, it is common to study fluctuation properties of the spectrum. In the 
present case this means that one disregards the unitary part of P„, t „ and passes to the following 
probability measure on C" (the density of states), 

(1.2) dn ffl ,„(A! A n ) = -^ |V„(Ai,...,A„)| 2 e- m ^U Q ^ ) dA n {A l ,...,A n ), 

/Limit 
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Here the partition function Z„, „ is given by 

(1.3) Z,„,„= f |y„(A 1/ ...,A„)| 2 e- mS H Q(A ' ) dA„(A 1/ ... / A„), 

JC" 

where we put 6A n (A\, ... , A„) = dA(Ai) ■ ■ ■ dA(A„); cL4(z) = d 2 z/7i is the suitably normalized area 
measure in the plane. (The integral | |1.3t converges when mjn > p~ l ; we always assume that this 
is the case.) 

Now fix a number x such that 

< T < p. 

We can think of the eigenvalues (A,)" as a system of point charges (electrons) confined to a 
plane, under the influence of the external magnetic potential Q |38]. In the limit when m —* oo, 
n/m — » t, the growth condition l jl.lt on Q is sufficient to force the point charges to condensate on 
a certain finite portion of the plane, called the "droplet", the details of which depends on Q and 
t. Thus the system of electrons, the Coulomb gas, lives in the vicinity of the droplet. Inside the 
droplet the repulsive behaviour of the point charges takes overhand and causes them to be very 
evenly spread out there. 

The droplet. We review some elements from weighted potential theory. Let A denote the nor- 
malized Laplacian, A = dd where d — \(d x — idy) and d = \{d x + id y ). Write 

(1.4) X = {AQ > 0}. 

Let SH T denote the set of subharmonic functions / : C — > M such that /(z) < Tlog + |z| 2 + 0(1) as 
z — > oo. The equilibrium potential Q T is defined as the envelope 

Q T (z) = sup {/(z); / e SH T , / < Q on C} . 
The droplet associated with the number t is the set 

(1.5) St = {Q = Q,|. 



Our assumptions then imply that Q T 6 SH T , Q T 6 C 1,1 ^), »S T is a compact set and Q T is harmonic 
outside 5 T . See e.g. |31| or [22J. In particular, since z i-> Tlog + (|z| /C) - C is a subharmonic 
minorant of Q for large enough C, it yields that 

(1.6) Q T (z) = Tlog + |z| 2 +0(l) on C. 

Let f be the convex set of all compactly supported Borel probability measures on C. The energy 
functional corresponding to t is given by 

There then exists a unique weighted equilibrium measure o T eP which minimizes the energy I T (a) 
over all a 6 !P. Explicitly, this measure is given by 

da T (z) = T- x AQ T (z) dA(z) = T- 1 AQ(z)l 5inx (z) dA(z). 

Cf. OH, 11221 . 

The problem of determining the details of the droplet are known under the names "Laplacian 
growth" or "quadrature domains". When Q is real-analytic in a nieghbourhood of the droplet, 
the boundary of the droplet is a finite union of analytic arcs with at most a finite number of 
singularities which can be either cusps pointing outwards from the droplet, or double-points, u) 
On the other hand, if Q is just C°°-smooth, the boundary will in general be quite complicated. 
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The correlation kernel. We state a couple of well-known facts concerning the measure Tl mn , 
l|1.2|l . For positive integers n with n < mp we let H m /H be the space of analytic polynomials of 
degree at most n - 1 with inner product {f,g) m Q = Lf( z )s( z ) e~"'Q (z) dA(z). We denote by K m/ „ 
the reproducing kernel for H m/n , i.e., 

n 
K min {z, W) = Y_ J <Pj(z)<Pj(w), 

where \<pj} n =1 is an orthonormal basis for H m „. 

It is well-known that H m ,n is given by a determinant 

(1.7) dn m ,„(Ai, . . . , A n ) = ±- det (K m , n (Ai, Aj) e -»(Q(A,) + Q(A,))/2\» ^ (A . . A „). 

More generally for k < n, the k-point marginal distribution IT^ H is the probability measure on 
<C k which is characterized by 

(1.8) f /(Ai A k ) dnf„ „(Ai A t ) = f /(And) A n(k)) dn ffl ,„(Ai A„), 

JC k JC" 

whenever / is a continuous bounded function depending only on k variables and n : {1, . . . , k] —* 
{1, . . . , n) is injective. Evidently Tl m ,n - njj, n . One then has that 

(1.9) dlt, n (h h)= { ^~ det (K m , n (A h Aj) e -™(Q(A,) + Q(A,))/2)^ =i ^ (Aj h) 

For proofs of the identities CLZ), CLD/ see e.g. (22, 122, or the argument in gl], §IV.7.2. 
The weighted kernel K m ^(z,zv) e~ m (Q( z ) + Q( IU ))/ 2 is known as the correlation kernel or Christoffel- 
Darboux kernel corresponding to the ensemble. 

Linear statistics. Let us now fix a function g e C;,(C) and form the random variable ("linear 
statistic") 

n 

trace,, g : C" ^ C , (Ay)? =1 » £ g(A ; ). 

Let E m „ denote expectation with respect to the measure n m „ on C" . Likewise, if k < n we let E k m n 
denote expectation with respect to the marginal distribution Yl k m n . Then by \\&\ , \\.9\ 

(1.10) E m , n (- trace,, g)=-f J fijgfrj)) = E^C^Aj)) = - f g(A x ) K m ,„(A!, A t ) er^^cLACAi). 

\n I n *r-f n Jc 

The asymptotics of the right hand side can be deduced from the following fact (see 11221 , cf . also 

0,GEG3) 



(l.ii) 



L 



- K m , n (A,A) e- nQW - t^AQAA) 
n 



dA(A) — » 0, as m — » oo, n/m — » t. 



Combining with l ll.lOt one obtains the following well-known result. 
Theorem 1.1. f[22J) Lef g e C fc (C). Then 

-E m „(trace„ ?) — > I g(A) dcr T (A), as m — > oo, n/m — > t. 
« ' Jc 

We now form the random variable ("fluctuation about the equilibrium"), 

fluct,, g = trace,, g - n J g dcx T . 
Jc 

The main problem considered in this paper is to determine the asymptotic distribution of 

fluct,, g as m — > co and « - m% —* when g is supported in the interior ("bulk") of S T n X. For this 
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purpose, we will use a result, due to Berman [7j, concerning the near-diagonal bulk asymptotics 
of the correlation kernel. 

Approximating Bergman kernels. For convenience we assume that Q be real-analytic in neigh- 
bourhood of the droplet. This is not a serious restriction, see Remark ll.5l and £17.11 (Moreover, 
the real analytic case is the most interesting one.) 

Let bo(z, w), b\{z, w) and ip(z, w) be the (unique) holomorphic functions defined in a neighbour- 
hood in C 2 of the set {(z,z);z 6 S T n X} such that b (z,z) = AQ(z), b\(z,z) = jAlog AQ(z), and 
i^(z, z) = Q(z) for all z e X. The first-order approximating Bergman kernel K} n (z, w) is defined by 

K] n {z,w) = (mb (z,w) + h(z,w)) e m ^ m \ 

for all z, w where it makes sense, viz. in a neighbourhood of the anti-diagonal {(z, z); z e «S T fl X}. 

Lemma 1.2. (|7J) Let Kbea compact subset ofS° D X, and fix z e K. T/zere £/ie« exz'sfs a numbers m , 
C and e > independent of zo swrfz that for all m>m§ holds 

| K m ,„(z, w) - K^(z, w) | g-«(QW+aw))/2 < cm-^ Z/ w 6 D ( Zo; e)/ n >mz-l. 
In particular, 



(1.12) 



K m ,„(z,z) e- mQ(z) - (mAQ(z) + ^Alog AQ(z)) 



< Cm , z e K, n> mx - 1. 



A proof of the result in the present form appears in [1], Theorem 2.8, using essentially the 
method of Berman [7\ and the approximate Berman projections constructed in |8| (compare also 
[12], [11]). Cf. \7\, §1.3 for a comparison with the line bundle setting. 

We remark that corresponding uniform estimates in Lemma 11.21 up to the boundary of the 
droplet, are false. 

Expectation of fluctuations. Using Lemma [l.21 we can easily prove the following result. 
Theorem 1.3. Suppose that g e Cg(S° n X). Then 

E m/n fluct,, g — » I g dv as m — > oo and n - mx — » 0, 
Jc 



where v is the signed measure 



dv(z) = ^Alog AQ(z) l 5Tnx (z) cM(z). 



Proof. By JTTl2t . 

E m ,„ (fluct,, g) = nE^,„g(Ai) - n I g(Ai) da T (Ai) = 

Jc 

= f (;«AQ(z) + iAlogAQ(z) + 6)( m - 1 ))^(z)dA(z)-«T- 1 f g(z) AQ(z) cL4(z) = 

Jsuppg \ -^ / Jsuppg 

= (m - m- 1 ) fg(z) AQ(z) cL4(z) + 1 f g(z) Alog AQ(z) dA(z) + 0{m~ l ). 
When m — » oo and ffz - «t _1 — > 0, the expression in the right hand side converges to JL gdv. D 
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Main result. Let V = (d/dx,d/dy) denote the usual gradient on C = ]R 2 . We have the following 
theorem. 

Theorem 1.4. Let g 6 C^ (S° n X). The random variable fluct,, g on the probability space (C", Yl mi „) 
converges in distribution when m — » oo and n - mx — > to a Gaussian variable with expectation e g and 
variance v 2 , given by 



h 



fgdv , v 2 g = ±J\Vg\ 2 dA. 



This theorem is the analog of a result due to Johansson [26 J, where the Hermitian case is 
considered. Following Johansson we note that in contrast to situation of the standard CLT, there 
is no 1/ V^-normalization of the fluctuations. The variance is thus very small compared to what 
it would be in the i.i.d. case. This means that there must be effective cancelations, caused by 
the repulsive behaviour of the eigenvalues. One can interpret Theorem 1 1.41 as the statement that 
the random distributions fluct„ converge to a Gaussian field on compact subsets of the bulk of the 
droplet. See 3731 

The formula for e g has already been shown. The rest of the paper is devoted to proving the 
other statements, viz. the formula for v 2 „ and the asymptotic normality of the variables fluct,, g 
when m — > oo and n — mz —* 0. A simple argument shows that it suffices to show these properties 
for real-valued functions g. In the following sections we will hence assume that g is real-valued. 

We will give two proofs of Theorem ll.41 one using the well-known cumulant method will be 
given in detail and another using an idea of Johansson |26| is sketched in 37.21 A comparison is 
found in Remark |7T1 

We here want to mention the parallel work by Berman [5J, who independently gave a different 
proof of a version Th. ll.4l valid in a more general situation involving several complex variables. 



Remark 1.5. We emphasize that in our first, cumulant-based proof of Th. I1.4l we assume that Q 
be real analytic in a neighbourhood of the droplet. The theorem is however true e.g. also for 
general Q : C — > IR U {+00} which satisfy l jl.1} and are finite and C°° except in a finite set where the 
value is +00. Since this type of potentials are sometimes useful, we will after the proof indicate 
the modifications needed to make it work in this generality. See 37.11 



The cumulant method. For a real-valued random variable A, the cumulants Ck(A), k > 1, are 
defined by 

00 ,k 

(1.13) logE( e M ) = £-C,(A), 

fc=i 

and A is Gaussian if and only if Ck(A) = for all k > 3. Moreover, C 2 (A) is the variance of A. 

It was observed by Marcinkiewicz that in order to prove asymptotic normality of a sequence of 
r.v.'s (i.e. convergence in distribution to a normal distribution), it suffices to prove convergence of 
all moments, or, equivalently convergence of the cumulants. Indeed convergence of the moments 
is somewhat stronger than asymptotic normality. 

We now fix a real-valued function g e CJf (S° n X) and write Cm,n,lcig) for the fc'th cumulant of 
trace,, g with respect to the measure n,„ „ . Following Rider- Virag 1 30 1, we can write the cumulants 
as integrals involving the cyclic product 

(1.14) R m , n ,k(M, . . ., A*) = K,„,„(Ai, A 2 ) X,„,„(A 2 , A 3 ) • • • K nhn {h, Ai) e -'»(Q(A 1 ) + ... + Q(A,))_ 
Namely, with 

d- 15 ) G ^ *>=E*=t- E rari^ 

7=1 * k 1 + ...+kj=k,ki,...,k j >\ I' 1=1 
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we have (]30|, cf. also ED, fl35), (36)) 

(1.16) C m ,n,kig) = \ G k (Ai, ..., At) R,„,n,k(M, . • . , A t ) dA, c (Ai, . . . , A k ). 

Note that if Gjt(Ai, . . . , At) ^ 0, then A, e supp g for some i. 

The representation jl.161 was used by Rider and Virag [30] in the case of the Ginibre potential 
Q = \z\ to prove the desired convergence of the cumulants. In another paper |29], the same 
authors applied the cumulant method to study some determinantal processes in the model 
Riemann surfaces, and they prove analogs of Th. ll.4l for a few other special (radial) potentials. 

The methods of |30|, |29| depends on the explicit form of the correlation kernel. In the present 
case, the explicit kernel is too complicated to be of much use. To circumvent this problem we will 
use the asymptotics in Lemma POl and also some off-diagonal damping results for the correlation 
kernels (cf. section |5). 

We want to emphasize that the result of [30] covers also the situation when g is not necessarily 
supported in the bulk (in Ginibre case), and this situation is not treated in Th. 11.41 (We shall have 
more to say about that case in general in 37.4l below.') 

The cumulant method is well-known and has been used earlier e.g. by Soshnikov |35| and 
Costin-Lebowitz Ifl6l to obtain results on asymptotic normality of fluctuations of linear statistics 
of eigenvalues from some classical compact groups. The method has also been used in the parallel 
work on linear statistics of zeros of Gaussian analytic functions initiated by Sodin-Tsirelson |34| 
and generalized by Shiffman-Zelditch |33|. A brief comparison of these results to those of the 
present paper is given in 37,81 



Other related work. It should also be noted that Th. 11.41 as well as the more general Th. 17.31 
below, follow from the well-known "physical" arguments due to Wiegmann et al. See e.g. the 
survey 1 38 1 and the references therein as well as 1 20 ] . 

Results related to fluctuations of eigenvalues of Hermitian matrices, are found in Johansson 
l26l and also 0, [4] and [21 j. A lot of work has been done concerning ensembles connected with 
the classical compact groups. See e.g. JT7), (2JQ, ||35l, (37J, (16J. 

Disposition and further results. Sections |2]-[6] comprise our cumulant-based proof of Theorem 
11.41 In our concluding remarks section, Sect. we state and prove further results. We summarize 
some of them here. In £ )7.2l we sketch an alternative proof of Th. ll.4b ased on a variational approach 
in the spirit of Johansson's paper [26 J. In 37.41 we state without proof the full plane version of 
Th. 11.41 The proof will appear in [2J. In 37.5| we prove universality under the natural scaling: if 

m = n then for a fixed zo 6 S^ n X, the rescaled point process ( Ay J i-» ( yfn(Aj - zq)\ _ converges 

to the Ginibre(oo) determinantal point process as n — > oo. In 37.6| we clarify the relation of our 
present results to the Berezin transform (which we studied in [1]); in particular we prove the 
"wave-function conjecture" ([22J) that |P„| e~ n ^dA converges to harmonic measure at oo with 
respect to C \ St, where P„ is the n:ih orthonormal polynomial corresponding to the weight e~ n Q 
and C = C U {co}. 

2. Further approximations and consequences of Taylor's formula. 

In this preparatory section, we discuss a variant of the near-diagonal bulk asymptotics for the 
correlation kernel and for the functions R n ,m,k (see (|1.14t ). especially for k = 2,3; such estimates 
are easily obtained by inserting the asymptotics in Lemma [l. 21 and they will be used in Sect. [6] 

In this and the following sections, we assume that Q is real analytic near the droplet, except 
when otherwise is specified. Recall that ip denotes the holomorphic extension of Q from the 
anti-diagonal, i.e. ip(z,z) = Q(z). 
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It is well-known and easy to show that ip is determined in a neighbourhood of a point at the 
anti-diagonal by the series 



_ . i jT j 

i ;=n J' 



for h and A: in a neighbourhood of 0. 

For clarity of the exposition, it is here worthwhile to explicitly write down the first few terms 
in the series for ip and Q 

ip (z + h, z + k) = 

= Q(z) + dQ(z) h + dQ(z) k+- (d 2 Q(z) h 2 + ^Q(z) F] + AQ(z) hk + "higher order terms", 

and 

Q(z + h) = Q(z) + dQ{z) h + dQ{z) h+ l - (d 2 Q(z) h 2 + ^Q(z) /?) + AQ(z) \h\ 2 + 0(\h\\ 



for small \h\. Using that ip(z, w) = ip(w, z), and that Q is real-analytic near the droplet, it is easy to 
prove uniformity of the (9-terms in z when z e S T . This means that there is £ > such that 

(2.1) |2Re^(z + / ; ,z)-Q(z)-Q(z + / ; ) + AQ(z)|^| 2 |<C|/z| 3 , z 6 <S t/ \h\ < e. 
We will in the following consider h such that \h\ < M5 m where M is fixed and 

5„, = log ml yfin. 
We then infer from | j2.1fr that there is a number C depending only M such that 

| 2m Re \p{z + h,z)- mQ{z) - mQ(z + h) + mAQ(z) \h\ 2 \ < Cmb\, z e S T , \h\ < M6 m , 

and mb\ = log 3 m/ yfm — » when m — » oo. 

Next recall the definition of the approximating kernel Kj,(z,w) = (mbo(z,w) + b\{z,w)')e m ^ z '' ui) 
(see Lemma [l.2t . We obviously have 

(2.2) \b (z + h, z) - AQ(z)| < C5 m and \h(z +h,z)\<C when z 6 5 T/ |/z| < M6 m/ 

for all large m with C depending only on K and M. It follows that 

(2.3) 

¥} m {z + ft, z) e -»<(Q(z+/0+Q(z))/2 _ m(A Q(z) + 6>(6 m )) e M^+Kz)-(Q(z) + Q(z + h))/2) / z e s ^ | fe | < M6m/ 

when m — > oo. Here the (9-term is uniform in z e »S T . Lemma ll.2l now implies the following 
estimate for the correlation kernel. 

Lemma 2.1. Fix a compact subset K c S° D X. Then for allz e K we have that 

\K m , n (z + h,z)\e- m ^ +h ^^ 2 = 

= m (AQ(z) + 0{b m )) e -'«AQWIftl 2 /2+0(iog 3 »/ V^) + 0{m~ l ), \h\ < Mb m , 

when m — > oo and n > m% - 1; the O-terms are uniform in zfor z e K. 

We will need a consequence concerning the functions R m ,n,k for k = 2 and k = 3. 

Lemma 2.2. Lef Kbea compact subset ofS° n X. Then for z e K we have 

(2.4) R ffl ,„, 2 (z, z + ft) = m 2 (AQ(z) 2 + 0(<5 m )) e - mAQ{z)m2+0 ^ 3 m/ ^ + 0(1), \h\ < Mb m , 
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and 

R m ,n^(z,z + hi,z + h 2 ) = 

(2.5) = m 3 (AQ(z) 3 + 0(5 m )) e ^Q^)(^-\h i \ 2 -\h 2 \ 2 ) + 0(iog i m/VJr,) + 

+ 0(l + m ( e -"'*Q(^\ 2 n + e -mAQ(z)|ft 2 | 2 /2J) / | fcl | / m < M5m/ 

when m — > oo and n > mi - 1; the O-terms are uniform for z e K. 

Proof. The estimate (|23) follows from Lemma|2l]sinceR m ,„,2(z / z+/i) = |K m , n (z,z + h)f e - m( - Q( - z)+Q< - z+K >\ 
To estimate R m ,n£(z,z + h\,z + hi) we first consider the approximation 

R\^Z,Z + h,Z + h 2 ) = K l m {z,Z + h)K^(z + h lr Z + h 2 )K l m {z + h 2 ,z) ^(QW+Q^HQfz+fe)^ 

obtained by replacing K„ hn by Kj n in the definition of R m ,„,3. 
In view of (12.21 we have for z e K and |fri| , |/z 2 | < Mb m that 

(2.6) R^z, z + ftj, z + h 2 ) = m 3 (AQ(z) 3 + 0{b m j) ^(#^M)+*(^^)+#+fe-2)-QW-QW)-Q(^)) ; 

where (9 is uniform in z 6 K. A simple calculation with the Taylor expansions for Q at z and ip at 
(z, z) now yields that 



i/>(z,z + hi) + i/>(z + hi,z + h 2 ) + i/'(z + h 2 ,z) - Q(z) - Q(z + h) - Q(z + h 2 ) = 
= AQ(z)(h 1 h 2 -\h 1 \ 2 -\h 2 \ 2 )+0(\h\i) r as fc^O, 

where we have put l^ = maxj|/ii| , \h 2 \}. Since the estimate is uniform for z e K, we may use | |2.6t 
to conclude that 

R^ 3 (z,z + hi,x + h 2 ) = m 3 (AQ(z) 3 + 0(6 m j) e "^Q(^h 2 -\h\ 2 -\h 2 \ 2 ) + 0(io g 3 m/^H / ^ < M6m/ 

when m — > oo and again the O-terms are uniform for z £ K. Combining with Lemma 12.11 and 
(l23t , and also using the estimate |X m ,„(z + /zi,z + /z 2 )|e" m(Q(2+, ' l)+Q(2+ ' !2))/2 < Cm for |/z| M < M5 m , 
n > mr -l,tn large (this follows from Lemma [L2t we readily obtain | |2.5| |. n 

3. The functions G k ; near-diagonal behaviour. 

In this section, we let g be any sufficiently smooth (sometimes real-valued) function on C (i.e. 
not necessarily supported in 5° D X). We then form the corresponding function G k by | |1.15t . Here 
k > 2 is fixed. 

We will now analyze the function G k in a neighbourhood of the diagonal 



where 



A k = {Al k eC k ;AeC], 



l t = (l # l / ... / l)6C*. 



Our results in this section state that G; c vanishes identically on Afc and that G k is harmonic at 
each point of A^. This depends on combinatorial identities of a type which where considered 
earlier in related contexts e.g. by Soshnikov [35J and Rider- Virag [30], [29 1. The following lemma 
is equivalent to [35J, eq. (1.14), p. 1356. 

Lemma 3.1. For any function g : C — » C and any k>2, it holds that G k = on A k . 

Proof. Evidently 



G k (Al k ) = g(A) k Y J { -^— Yj 

7=1 ' h+...+kj=k,k 



-* ki\---h\ 

;h,...,kj>\ I 
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The last sum is the number of partitions of k distinguishable elements into / distinguishable, 
nonempty subsets. Thus (e.g. ED, Th. 9.1, p. 340) 



L 



k\ 



where 



i—i h, I . . . fc-l 

k l +-+kj=kjc u ...,k } >l ! ' 



flS(k,j), 



is the Stirling number of the second kind. Evidently S(k,0) = for k > 1. Moreover, the 
well-known recurrence relation for those Stirling numbers (see e.g. [14], Th. 8.9, (8.32)), gives 



/,--! 



(-l)H 



S(k - 1, 0) = £(-l)'r!S(*, r + 1) = £ ^J 

r=0 ;=1 ' 



-jlS(k,j). 



The lemma follows, since S(k - 1, 0) = when k > 2. 
Note that the lemma is equivalent to that 



(3.1) 



(-l)i- 1 






fcil-.-fc! 



= 0, jfc = 2,3,.... 



it l +»-+j^=iyti^..^i " 1 ' ' ' 
We note the following simple, but rather useful consequence of Lemma |3T 
Lemma 3.2. Let g e C : (C ^ C) and k>2. Then for allAeC holds 

k k 



!=1 



T 0A)(Ai, • • . , A,) = y (5,-G t )(Ai, . . . , Ait) 



/=! 



= 0. 



=At=A 



Proof. By Lemma [3.1 1 we have that G)t(Alfc) = 0, whence 

k 



= 4rG k {M k ) = YjdiGk){Mk). 

The statement about 5 is analogous. □ 

We now turn to a more nontrivial fact. Let us denote by 

A k = d\d\ + . . . + d k dk, 

the Laplacian on C k . 

In the next lemma, we calculate A k G k at every point of the diagonal A^ when k > 2. When 
k > 3, we shall see that A k G k vanishes on the diagonal, which means that G k is nearly harmonic 
close to the diagonal. 

Lemma 3.3. Let g 6 C 2 (C -> R) and k>2. Then for all A e C we have 



(A 2 G 2 )(A 1 ,A 2 ) 



Ai=A 2 =A 



Vg(A) 12, 



and 



(A k G k )(A,,...,A k ) 



= 0, fc = 3,4,... 



A 1= ..=A k =A 
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Proof. Fix a number k > 2. Let 1 < j < k, and let ki,...,h be positive integers such that 
ki + . . . + kj = k. Since, for 1 < r < /', 

~o I i \ i _ i 

^pll^' =*>(*,-!)• I] g^i) k '-g^r) k '- 2 -dg(A r ydg(A r )+k r - Yl giAtt-giArf-'-AgiAr), 



.1=1 



l=U*r 



\=l,l±r 



we get (with l k = (1, . . ., 1) 6 C*) 



(A fc 



(-l)i- 1 



(3.2) 



y=l ; &i+...+Ji: / =jyci,...,fc;>l ' 

2 i 
g(A)*" 2 <?g(A)| £ fc(*r - 1) + ^A)*" 1 Ag(A) £ fc r 



/ ^ 



Since k\ + . . . +kj = k, the right hand side in | |3.2| | simplifies to 



2 A (-l)H 



(3.3) 



^wiYFi- E 

/=1 I ki+...+h=k,k 



fc!(fcl(fcl-l) + ... + fc ; (fc;-l)) 



&!+...+£,=£, £!,...,£,■>! 



h\---kj\ 



;=1 ' ki + ...+kj=k,k 



k-kl 



-> k l \---k j \ 

k,h,-,kj>l ' 

Here the last double sum is zero, by fl3.lt , and (|3.3t simplifies to 

I 2 V H) 7 " 1 X- 1 fc! ( fc i( fc i -!) + •■• + */(*/ - 1)) 



(3.4) 



sw k - 2 d g (4Z^f- E 

;=1 ' h+...+k.=k,k 



;-'W»"/'V = 



hl---kjl 



In order to finish the proof we must thus show that S2 - 2 and Sjt = for all k > 3 where S,t 
denotes the sum 



(3.5) 



(-l)H 



E 1 ^ E 

;'=1 ' h+...+kj=k,k 



ki+...+kj=k,ki,...,kj>l 



fcl!---fcy! 



The case fc = 2 is trivial, so we assume that k > 3. To this end, we shall consider exponential 
generating functions of the form 



(3.6) 



H 



• ~ fci=i '■ ■ " 



i=i 



k; = l 



The relevance of this generating function is seen when we expand the product as a power series 
in t, 



( 



1 1 



H / (f;xi ^ )= E E 

k=l \ki+...+kj=k,h kj>l I 1 



Considering the Xf.S as real variables and denoting 



A = h . . . H , 



dx 2 



the Laplacian on W, we thus obtain 



(3.7) AfH 7 (f;l 1) = £ £ 

\k- l +...+ki=k,ki,...,kj>l 



dx 2 , 



k\(ki(h-i) + ...+kjgcj-i)) 



k=l 



h\---kj\ 



M 
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On the other hand, differentiating the product in (13. 6) and evaluating at X\ 



1 yields 



(3.8) AfHj(f;l,...,l) = /W-l) / ~ 1 , 

Differentiating J3.7) k times with respect to t and evaluating at t = 0, we obtain the result that 

ttfafa -1) + ... + *;(*;•-!)) d* 



E 

fci+...+fcy=fcfci,_.,fc,'>i 



h\---kf. 



In view of (13. 5t , this implies that 

( k 

(3.9) 



s - d * 
s *"d7 



=S(^-r) 



^M-4)) 



f=0 



7=1 
But since 1 - e t = —(t + t 2 /2\ + f 3 /3! + ...), it is seen that the coefficients fl; in the expansion 

CXI 



1=0 

must vanish whenever I + 2 and I < k + 2. In particular, if, as we have assumed, k is at least 3, 
then we have a k = 0, which by | |3.9t implies that S k = 0. The proof is finished. □ 

In addition to the Laplacian (AjtGfc)(Aljt), we will also need to consider functions of the form 

(3.10) Z k (A) = Yjd^G^iA.%), k > 2. 

The following lemma is now easy to prove. 

_ 2 

Lemma 3.4. We have that Z%{?C) - - <?g(A) while Zj is pure imaginary when k>3. 

Proof. Again the case k = 2 is trivial because G2(Ai, A2) = g(Ai) 2 - g(Ai)^(A2). When k > 3 we may 
use lemmas |3~T1 and 1331 to calculate 



= A A {G k (Al k )} = (A k G k )(Al k ) + Yj.didjG k )(Al k ) = 2ReZ k (A), 

i*i 

which shows that Z k is pure imaginary □ 

4. An expansion formula for the cumulants 

During this section, we keep a real valued function g 6 C™(S° n X) fixed. We will reduce the 
proof of Th. 1 1.41 to the proof of another statement (Th. l4.4l below), which turns out to be easier to 
handle, and which we prove in the Sect. [6j after a discussion of some basic estimates for K m r „ in 
Sect.0 

To get started, note that an expression for the cumulant C mj n,k{g) wa s given above in eq. dl - 16b - 
It will be important to note that ( 11.161 1 and the reproducing property of X m „ shows that we may 
also represent the cumulant C m ,n,k{g) as an integral over C k+1 , 



(4.1) 



C m , n ,k(g) = Gjt(Ai, . . . , A k ) R mi n lk +\{A, A\,..., A k ) dA fc+1 (A, Ai, . . . , A k ), 



where G k and R„ h „, k +i are given by dl,15t and l|1.14t respectively. Indeed, this simple trick of 
introducing an extra parameter A into the integral will turn out to be of fundamental importance 
for our proof. 

In the foregoing section, we were able to give a good description of Gjt(Ai, . . . , A k ) for points 
near the diagonal A\ — . . . — A k — A. For such points it is natural to write hj — A,- — A (where the 
\hi\ are small) and to work in the coordinate system (A, h\, . . . , h k ). Indeed, this coordinate system 
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is advantageous for all our purposes. Note that the volume element is invariant with respect to 
this change of coordinates, 

dA k+1 (A,Ai,...,A k ) = dA k+1 (A,h,---,h), 
and that the reproducing property of K m r „ is reflected by the fact that 

u(A) = f u(h) K m , n (A, A + h) e - mQ{A+h) dA(h), u e H m , n . 
Jc 



We thus get that with h = (hi, . . . , h k ) and l k = (1, . . . , 1), we can write (|4.1|l as 

(4.2) C m , n , k {g)= \ G k (Al k + h)R m , hk+1 (A,Al k + h)dA k+1 (A,h). 

Jc k+i 

We now fix A e C and use Taylor's formula applied to the function 

C k -> K : )ih G k (Al k + h). 
Since G k (Al k ) = by Lemma [3. 1[ the Taylor series at h = can be written 

(4.3) G k (Al k + h)~Yj T M> h )' 
where, in the multi-index notation, 

[d^ G k ){Al k )—. 
|«+j5|=; 
Note that if A £ supp g, then G k vanishes identically in a neighbourhood of Al k , and so Tj(A, h) = 
for all h e C*. Thus the right hand side in J4.31 is identically zero when A £ supp g. 

Let us write |/z|oo = max{|/ii| , . . . , \h k \}. It will turn out to be sufficient to consider Taylor series 
of degree up to two. We thus put 

(4.4) G k (Al k + h) = T 1 (A,h) + T 2 (A,h) + r(A,h), where r(A,h) = 0(\h\i) as 7z->0. 



The idea is now to replace G k (Al k + h) by the right hand side in (|4.4) in the integral l|4.2) . To 
simplify matters, we first have the following lemma. 

Lemma 4.1. For allk>2 holds 

Ti{A,h) R nhnM1 (A,Al k + h) dA k+1 (A,h) = 0. 



f 

Jo- 



Proof. First note that 



(4.5) Ti(A, fe) = 2 Re £](5A)(Alt)fc,-. 



Integrating termwise in (14. 5t with respect to the measure R„ h „ rk +i(A, Al k +h)dA k+ i(A, h) and observ- 
ing that the terms on the right hand side of J4.51 depends only on two variables, the reproducing 
property of K,„ r „ shows that, for i = 1, . . . , k, 

\ (d t G k )(Al k ) R m , n , k+1 (A,Al k + h) hi dA k+1 (A,h) = \ (d,G k )(Al k ) R m ,„, 2 (A,A + h) h x dA 2 (A,/n), 
Jc ,c+1 Jc 2 

and so we can replace the integral in (14. 5t by an integral over C 2 (since R m , n ,2 is real-valued): 



f T 1 (A,/z)R„,, a+1 (A, Al t + ft)d4t+i(A,ft) = 2 Re f f'(aA)(Al*) 

Jc t+1 Jc 2 (^ 

The last integral vanishes by Lemma [3. 21 



R m ,n,2{A,A + hi)h l dA 2 {A,h l )- 
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We have shown now shown that 



C m ,„, k (g) = \ (T 2 (A,h) + r(A,h)) R m , nM i(A, Al k + h) dA k+1 (A,h). 

Jc k+l 



To simplify this expression further, we will first look more closely at 



(d"J G k )(Al k )— , 

|«+/3|=2 



which we write in the form 



w=i 



1 k 1 k k _ 

T 2 (A,h) = - J^OidjGkXAWhjhj + 2 J^idid^GkiAWkhj + JjdidjG k ){Al k )h 
i,7=l ',7=1 

k 

= ReJ^(d*G k )(Al k )h* + ReJ^(d i d j G k )(Al k )h i h j + 

i=\ i*j 

k 

+ J^(dJiG k )(Al k ) \h,\ 2 +2ReJ^(did j G k )(Al k )h i h j . 



,hj 



;=i 



'<] 



Using the reproducing property of X m „, it yields (note that R m ,n,k is not real-valued if k > 3) 



r t 2 (a, 



^)R m ,„,jt+i(A,Aljt + h) dA k+ i(A,h) 



( 



= f Re V(^;G k )(A1*)/»i/» 2 

Jc3 It? 

+ Re | I r,(^Q)(Al,c) 



\ 



R m ,«,3(A, A+h,A + h 2 ) dA 3 (A, hi, h 2 )+ 



1(5 

r ( 

2 I Re ^(^ y Gt)(Aljt)/»iS2 
V ((^,G,)(Al,))|/7i| 2 R m , n/ 



hlR m , n ,2(A,A + hi) dA 2 (A,hi)+ 
\ 



Rm,«,3(A, A + hi,A + h 2 ) dA 3 (A, /ii, h 2 )+ 
, 2 (A,A + hi)dA 2 (A,hi). 



Let us now introduce some notation. Recall that 



(AjtQtXAl*) = ^(^GtXAljt) and Z*(A) = J^(didjG k )(Al k ), A e C. 

;=i /<; 
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Definition 4.2. Let us put 

Rm.nM^ A + h, A + fc 2 ) dA 3 (A, fti, h 2 ), 



A m/ „(k) = f Re T'(^yGt)(Al*)ftift 2 

Jc3 It/ 



IS ( 



B m , n (k) = Re > .(^GtXAli) 



^R m ,„, 2 (A,A + fci)cL4 2 (A,/ji), 



C m , n (fc) = 2 f Re (z k (A)hh 2 ) R m ,„, 3 (A, A + h lt A + fe 2 ) dA 3 (A, h x , h 2 ), 

D m , n (k) = \ (A k G k )(Al k ) \h\ 2 R m ,„, 2 (A, A + h x ) dA 2 (A, hi), and, 
Jc 2 

E m ,„(k) = I r(A, h)R m , nik+1 (A, Alt + h) dA t+ i(A, ft). 

Jc t+1 

Our preceding efforts in this section are then summed up by the following formula. 
Lemma 4.3. For all m, n, k and all g e C^°(C) we have 
(4.6) C,„,n,k(g) = A m>n {k) + B m ,„(k) + C m ,„(k) + D m , n (k) + E m , n (k). 

The rest of this paper will be devoted to a proof the following theorem. 

Theorem 4.4. Suppose that g e C^(S° n X). Then for allk>1 the numbers A m/ „(k), B m/ „(k), and E„ h „(k) 
converge to as m — > oo and n - nn — » 0. Moreover we have that 

lim D m , n (fc) Ji/clW)| 2 cM(A) if k = 2, 

m-»oo,n-mT-»0 |0 if k > 3, 

and 

lim C ffl ,„ W = {-^clW)| 2 dA(A) Sf «: = 2, 

m— *oo,H-mT— >0 I Q if k > 3. 



It should be noted that Th. l4.4l implies Th. ll.4l (Convergence of the cumulants of fluct,, g to the 
cumulants of N \e g , vVj is equivalent to convergence of the moments which implies convergence 
in distribution.) 

In order to verify Th. 14.41 we will first need to look more closely at the behaviour of the 
function (A, h) i-> Gj^Aljc + ft)R m ,n,/c+i(A, A\ k + h) in the next section. We shall see that this function 
becomes negligible when h is "large" in the sense that |ft,| > M k log mj y[m for some i, where M k is 
a sufficiently large number independent of m and n as long as supp g C S° fl X and \n - m%\ < 1. 
This will imply that we can approximate the integrals defining the numbers A„ hn (k),. . . , E mi „(k) 
by integrals over a small neighbourhood of the diagonal in C k+1 . 

5. Off-diagonal damping 
Fix a number k>2. Throughout this section, it will be convenient to denote 

Ao = A k+ i = A, 
so that we can write 



R m ,n,k+\{A, . . ., Afc) = I I -K„,,„(A„ A,+i) e 



m(Q(Ai)+Q(A i+1 ))/2 

i=0 



We will frequently without further mention apply this convention in the sequel. We will need 
two lemmas. 



X 
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Lemma 5.1. ([7]) There is a number C such that for all z,w eC and all m, n with n < mx + 1 holds: 

\K mn (z,w)\ 2 e -"'(Q(z)+Q(a')) < Q m 2 e -m(Q(z)-Q z (z)) £ -m(Q(w)-Q 1 (w)) _ 

Proof. See Q or |T|, Prop. 3.6. n 

Lemma 5.2. (HI) Let Kbe a compact subset ofS° n X and d = dist (X; C \ (S T D X)). TTzere f/zen exist 
positive numbers C and e depending only on d such that for all z e K,h e C and a// m,n>\ such that 
\n - mx| < 1 ZioZds: 

Proof See Q, Th. 8.3, cf. also 0. n 

It follows from Lemma |5H that 

(5.1) \R m ,„, k+ i{A, Ai, . . . , A k )\ < Cm k+1 r»«2W-Q.W) ^(OC^-S^)) . . . g^OW-atA*)), 

when n < mx + 1. By the growth assumption l jl.lt , using that x < p and eq. (II. 6t , we conclude 
that there exists positive numbers C, C and 6 such that 
(5.2) 
|Km,n,ic+i| < C'm k+1 (max{|A| 2 ,... , \A k \ 2 }) ""' when n < mx + 1 and max{|A| 2 ,. . .,|A fc | 2 } > C. 

Thus if Dc(0) denotes the polydisc {(A, . . . , A k ); max{|A| 2 , . . . , |Ajt| 2 } < c}, we have for any NeE 
(|A| 2 + . . . + |A,t| 2 ) \ R mA/k+1 (A, . . . , A k ) | dA k+1 (A, . . . , Ak) -» 0, as m->oo, n<mz + l, 

C t+1 \D c (0) 

when C is large enough. We shall now show that much more is true. We first have the following 
lemma. In the proofs we conform to previous notation and write 

b„, = log m/ Vm. 

We also put 

d = dist (supp g; C \ (S T n X)) , 
and 

(5.3) K = {z 6 C; dist (z; C \ (S T n X)) > d/2} . 

We also remind the reader of the convention that A k+ \ = Ao = A. 

Lemma 5.3. TTiere exists positive numbers M, a and mo depending only on k and d such that if A; 6 X 
and Ay - Ay+i > M5 m for some index j 6 {0,. . .,fc}, then for all m > mo 

| £m,n,jfc+i(Ao, Ai, . . . , Ajfc) | < Cm" tt , |n - mx| < 1, 

o;/iere C depends only on d. 

Proof. In view of Lemma I5T21 the hypothesis yields that 

| K m , n (Aj, Aj +1 ) | e -™(Q(A / )+Q(A ;+1 ))/2 < Cm e ,-eVSmta{i/Z,|A r A /rt |} / |„ _ m% \ < ^ 

with numbers C and e depending only on d, and Ay - Ay+i > MS m . Choosing mo large enough 
that M6 m < djl for m > m it yields that 

(5.4) | K m , n (Ap A /+1 ) | e -'»(Q(A,)+Q(A /+ i))/2 < Cm e -eyfcM6 m _ Cm 1 " 6 * 1 , |n - mx| < 1, 
when m > mo. On the other hand, if n < mx + 1, Lemma [5.1| yields that 

(5.5) I K mtn (A u A J+ i) I e -™(Q(A ( )+Q(A; + i))/2 < Cffl/ ; _ 0/ . . . , fc. 
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Now i5A\ and (15. 5t implies 

k 

(5.6) | R m , nM1 (A A k ) | = J] | K m>n {A lr A /+1 ) | e -<"(Q(A,)+Q(A, +I ))/2 < C m M - M , 

1=0 

when m> niQ and \n - mz\ < 1. It now suffices to choose M large enough that 

eM - k - 1 > 0, 
and then put a = eM -k—1. a 

We henceforth let M denote a fixed large number with the properties provided by Lemma IB~3l 
Let us also put 

Ug(A) = dist (A; supp g) , A e C, 

U* S (A , . . . , A k ) = max { U g (Ai); i = 0,...,k], 
and 



y m/t = (^(A ,...,Afc)>MfcS m }. 



Lemma 5.4. The function 

(5.7) (A , Ai, . . . , Ak) h> Gjt(Ai, . . . , A*) -R m ,„,ic+i(A , Ai, . . . , A*), 
converges to zero uniformly on the set v/ m/ / c asm-^oo and |n - mz| < 1. 

Proof. Since G^ is bounded, it suffices to prove that R m , n ,k+i converges to zero uniformly on the 
set 

V 'm,k = V m,k n SU PP G/c- 

Here we regard G/ c as a function of the variables Aq, . . . , Ak, which is independent of the parameter 
Ao- It is then clear that 

supp Gk c {(Ao, . . . , Ak); Ao e C, and A, e supp g for some i = 1, . . . , k) . 

Thus if (A , . . . , Ak) e V' m , , then there exists an index ? 6 {1, . . . , k] such that A, e supp g. Since the 
function R m/ „,fc+i(Ao, . . . , A/ c ) is invariant under the cyclic permutation i — > 1 •— >■ . . . i — > A: )— > of the 
indices, we can iv.l.o.g. assume that i = 1. Then, since U g (A\) = and ITJAi,. . .,Ak+i) > Mk5,„, 
there must exist an integer j e {1, . . . , k) such that |A; - A;+i| < Mb m for all indices / with 1 < I < j 
and Aj - Ay+i > M5 m . It then follows from the triangle inequality that 

(5.8) Ug(Aj) < \Aj - Ai| < Mfc6 m . 
If m is large enough that 

(5.9) Mkb m < All, 

then | |5.8t implies that Ay belongs to the compact set K (see (15.31 ) and [Ay - Ay+i| > M6 m . Hence 
Lemma |5~3l yields that 

|R m ,»a+i(A ,...,A fc ) I < Cm~ a 
for large m when \n - mz| < 1, where a > 0. This proves that R m , n ,k+i converges uniformly to on 

V ,. □ 

Let us now put 



N(A , . . ., A k ) = max |A,- - A i+ i| . 

0</<J: 



We shall next prove that the function GkR m ,n,k+\ is uniformly small on the set 

W m , k ■- {(Ao, . . . , A k ); U g (A , ...,A k )> Mkb m or N(A , ...,A k )>M6 m }, 
where M = M(k, d) is a number provided by Lemma [5.41 
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Lemma 5.5. The function 

(5.10) (A , Ai, . . ., Ajt) i-» G;t(Ai, . . . , At) K m ^+i(Ao, Ai, . . ., At) 
converges to zero uniformly on W m ^ asm ^ oo and \n - mx| < 1. 

Proof. By Lemma 15.41 we know that the function J5.101 converges to zero uniformly on the set 
{U* > Mkb m }. It thus suffices to show uniform convergence on the set 

W ' m ,k = {u* g (A ,...,A k )< Mkb m and N(A A k ) > M8 m } . 

Now note that if m is large enough that Mkb m < d/2, we will have 

with K as in (|5.3|) . Hence if (Ao, . . ., A k ) e W . , we will have that A, e K and |A, - A !+ i| > Mb m 
for some z. It then follows from Lemma [5.31 that R m ,«,fc+i(Ao/ ■ ■ ■ / Ajt) < Cm~" when \n - m%\ < 1, 
where a > 0. It follows that R mA/ k+i ~ * uniformly on W , , and the lemma follows. □ 

It is now advantageous to pass to the coordinate system (A, h) where A = Ao and hi = A/ - A for 
i — 1, . . . , k. Let us put 

|/tL=max{|fc,-|; 1 < i < *}, 
and 

(5.11) Y m ,, = ((A, ft) e C t+1 ; LZ g (A) < Mic6 m/ |fc| M < Mfc6 m ) . 

As we shall see, everything interesting goes on in the set Y m ^ when m is large and \n - mz\ < 1. 
Lemma 5.6. The function 

(A, ft) h-» G t (Alt + ft) R m ,„,i+i(A, AI* + ft) 
converges to zero uniformly on the complement ofY m ,k «sm->oo and |n - mx| < 1. 



Proof. In view of Lemma [5.51 it suffices to prove that if (A, ft) is in the complement of Y m ^, then 
(A, Ai, . . . , Ak) belongs to W m ^, where A, = A + ft;. But if (A, ft) g Y,„ ^ then either U g (A) > Mkb m , or 
| A - A;| > Mfc6 m for some i = l,...,k. But the latter inequality can only hold if | Ay - A /+ i | > Mb m 
for some ;', whence N(A, A\,..., Ajt) > Mb m . Thus, in either case, we have (Ao, . . . , Afc) 6 W m ,j; and 
the lemma follows. □ 

The following result sums up our efforts in this section, and is what is needed to prove the 
asymptotic behaviour of the cumulants in the next section. 

Lemma 5.7. We have that 

' G k (Al k + ft) R m , n ,k+i{A, A\ k + h)) | dA fc+ i(A,ft) -> 0, 



f 

Jc 



asm — > oo and \n - m%\ < 1. 

Proof. It follows from | |5.2| | that the integrals 

Im = I G k (Al k + ft) R m , n , k+ i(A, A\ k + ft) dA)t +1 (A 7 ft) 

Jc t+1 \D c (0) 

converge absolutely for large enough m and C if ft < mi + 1, and J„, ->0asm->oo and ft < mi + 1. 
The statement now follows from Lemma [5.6l P 
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Remark 5.8. Suppose that P(A,h) is a measurable function on C k+1 such that (i) P(A,h) = when 

A £ suppg and (ii) |P(A, fr)| < CM + \h\ I for some constants C and N. (We write \h\ for the £ 

norm on C*, so that WL < W 2 < k \h\i-) 
As above, we can then conclude that 

(5.12) J P(A,h)R mArk+1 (A,A1k + h)dA k+1 (A,h) ^0, as m -> oo, \n-mz\ < 1. 

Indeed, (|5.12t follows from Lemma 15.31 if we use also the estimate j6.1|l to estimate the part of 
integral over \h\ > C for C large enough. The details of a proof parallel our proof of Lemma [5.71 
but are simpler in the present case, since U g (A) = when P + 0. 

6. Conclusion of the proof of Theorem 14.41 



In this section, we prove Th. 14.41 As we have observed earlier, this theorem implies Th. 11.41 
and thus the story ends with this section. 

Our proof will be accomplished by estimating the various terms in the identity 

C m ,n,k{g) = A mtn (k) + B m>n (k) + C„ h „(k) + D m , n (k) + E„,,„(fc), 

see (14. 6) . We start by considering the "error-term" 

E m ,„(k) = r(A, h) R„ hnM1 (A, Al k + h) dA k+1 (A, h), 

where r(A,h) is the remainder term of order 3 from Taylor's formula applied to the function 
h i-> G k (Al k +h) at h - 0, see (pL4) . We have that r(A,h) = G k (Al k + h) - P 2 (A,h) where P 2 is a 
polynomial of degree 2 in h with the property that Pz(A, h) = when A £ supp(g). It follows from 
Remark [5T8] that, when m —* oo and \n - mi\ < 1, (with Y„ hk as in l|5.11t ) 

(6.1) f P 2 (A, h) R m ,„ ik+1 {A, A\ k + h) dA k+1 (A, h) -» 0. 

Jc t+1 \y„,, t 

Using I l6.lt and Lemma 15*71 we conclude that 

r(A,/z) R m/nM1 (\,Al k + h) dA k+1 (A,h) -* 0, 



( 

Jc*- 



when m — > oo and |n — mz\ < 1. In order to estimate the integral over Y mk , we first introduce 
some notation. 

For a measurable subset O c C N , let us denote the (suitably normalized) complex N-dimensional 
volume of U by Volw(O) = J Q dA N (Ai, ...,A N )- When N = 1 we write Area(Q) in stead of Voli(Q). 

For large m, the set Y mk is contained in the set 

whence 

Vol k+1 (Y m>k ) < Area(5 T ) (Mk5 m ) 2k = CS~?, 
with C a number depending on k, M and t. Furthermore, l|5.Hl yields that 

| "W+i(A, Al/c + fe)) | < Cm fc+1 , n < tot + 1, 

for all A and ft. Now, since |r(A, ft)| < C \h\ 3 < Cb 3 m when |/i| < Mk5,„, it yields 

r(A,ft) R m ,„, k+ i(A,Ah + "i)| dA, +1 (A,/i) < C6^, m t+1 Vol, +1 (Y m ,i.) = 



1 



= Cm fc+1 6^ +3 = Clog 2 ^ 3 ml Mm. 

Hence also the integral over Y Mik converges to when m — > oo and \n - m%\ < 1. We have shown 
that E„ h „(k) — » as m — > oo and |n - tot| < 1. 
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We next consider the term 

D nhn {k) = f (AfcQtXAljt) \h\ 2 R m ,„, 2 (A,A + hi) dA 2 (A,h)- 
Jc 2 

In view of Lemma [3. 31 we plainly have 

D m/n (k) =0 if k > 3. 
It thus remains to consider the case k = 2. In this case, Lemma [33] implies 

D„,,„(2) = i f |Vg(A)| 2 |ft| 2 R m , n , 2 (A, A + ft) dA 2 (A,ft). 
^ Jc 2 

It is clear from Remark [5781 that 

(6.2) f |Vg(A)| 2 |ft| 2 R„,,„, 2 (A, A + ft) dA 2 (A,ft) -> 0, 

J|/i|>2M5,„ 

asnuoo and \n - ott| < 1. To estimate the integral over {|ft| < 2Mb m ) we apply the asymptotics 
for R,n f n,i from eq. (|2.4t (with the compact set X replaced by suppg). It yields that there are 
numbers v„, converging to 1 when m —* oo such that 

f |V£(A)| 2 |ft| 2 R„,,„, 2 (A,A + ft) dA 2 (A,ft) = 

, . J|/i|<2M6„, 

(6.3) 

= ^ m m 2 |Vg(A)| |ft| 2 (AQ(A) 2 + 0(6 m j) e- mA ^ A W 2 dA 2 (A,h) + o(l), 

when m — > oo and n > mx—1. Now, for a fixed A e supp g, the change of variables E, - -\/mAQ(A)h 
shows that 

f (mAQ(A)) 2 |ft| V mAQ(A)l ' i|2 dA(ft) = f \E,\ 2 e~^ dA{E) -» 1, 

J|ft|<2MS„, J|£|<2Mlogm 

as m — » oo. Hence it follows from (16.21 and ||6.3t that 



D m ,„(2)-»1 J|v^(A)| 2 dA(A), 



/C 
as m — » oo and |« - mz| < 1. 

The complete asymptotics for D,„ rn (k) has now been settled, and we turn to the term 

B m , n (k) = Re J S(A) ft 2 R„,,„, 2 (A, A + ft) dA 2 (A,h), 

Jc 2 

where we have put 

i- 
S(A) = £$Gt)(Al*). 

1=1 

Note that supp S c supp g. Using Remark |5lfl we obtain (as before) that 

\ S(A) h 2 R m ,„, 2 (A, A + h) dA 2 (A, ft) -> 0, 

J\h\>1Mb m 

asm —> co and \n - mz\ < 1. When \h\ < 2M5 m we again use the asymptotics in l|2.4|l , which yields 
that there are numbers v m converging to 1 as m — » oo such that 



(6.4) J| ^ 2M6 



f S(A) ft 2 R,„,„ /2 (A, A + ft) dA 2 {A, ft) = 

J\h\< 2 M6„, 

f S(A) ft 2 (AQ(A) 2 + 0{b m )) e -' nA Q^ 2 dA 2 (A, ft) + o(l). 

J\h\<ZM&„ 



v m m 2 
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Now, using that, for a fixed A 6 supp g, 

( (mAQ(A)) 2 h 2 e-' nAQ{m2 dA(h)= f ^e^dA^) = 0, 

J\h\<2M6„, J|<J|<2Mlogm 

we infer that B„ hn {k) — » for all fc > 2 as m —> oo and |« - mx| < 1. 

There remains to estimate the terms A mA {k) and C„ hn {k). These terms are a little more compli- 
cated than the previous ones since they are defined as integrals over C 3 and not over C 2 . We first 
turn to the term A mn {k) which we now write in the form 

A m , n (k) = \[ {TWhh 2 + Wjhfo) R m ,„, 3 (A, A + h, A + h 2 ) dA 3 (A, h, h 2 ), 

where we have put 

T(A) = Jjd,djG k ){Al k ). 

It is clear that supp T C supp g. Furthermore, using Remark 15.81 we see as before that, with 

h = (h 1/ h 2 ) and \h\ x = max{\h\,\h 2 \}, 

f Re (T(A)hh 2 ) R m ,„, 3 (A, A + h lf A + h 2 ) dA 3 (A, h, h 2 ) -> 0, 

asm^oo and \n - mz\ < 1. When [h]^ < 3M5,„, insert the asymptotics for R„,,„ /3 provided by eq. 
(12.51) . It shows that there are numbers v m converging to 1 as m — » oo such that 

f T(A) h h 2 R m ,„t(A, A + h, A + h 2 ) dA 3 (A, h) = 

J\h\ oc <3M5 m 

= m 3 v m f T(A) h t h 2 (AQ(A) 3 + 0{b m )) e mA Q^ h ^-^ 2 -^^ dA 3 (A, h) + o(l). 

J\h\„<3M5„, 

Now fix A e supp g and put E,\ = -\jmAQ(A)hi and E, 2 = -\jmAQ(A)h 2 . We then have that 
f T(A) (AQ(A) 3 + 0(d m )) h h 2 e »'AQ(A)C>iS 2 -|ft 1 | 2 -i; i2 | 2 ) dA ^ _ 

= T(A) f (1 + 0(6,„)) l x i, 2 e^2-\Zi\ 2 -\Z2\ 2 dA ^y 

J|4L<3Mlogm 

Thus when we can prove that / = and /' = where 

(6.5) /= f £ifce 4llH4l|2 - |42|2 dA 2 (£ 1/ £ 2 ) and /' = f & | 2 «*&-l€ifM&f< dA 2 (&,&) 
Jc 2 Jc 2 

we will obtain the result that A„ hn (k) — » as m — » oo and |n - mr| < 1 for all fc > 2. 

The argument for /' is similar so we settle for proving that / = 0. To this end, we write the 
integral in polar coordinates: 

J= — I I I(r,p)drdp, 

71 Jo Jo 

where 

Jr*2n pin 
I (rp) 2 e^ ^ e ^ iia -' p) -r 2 -p 2 d( p d0 . 
o Jo 
Performing the change of variables 9 = 6 + n/2 and <p = cp + n/2, the latter integral transforms to 

Jr.2n pin 
{rpf e i{n+9+ ^ e r P ^-r 2 -p 2 dd d(p = _ J(n p) 
o Jo 
Hence I(r, p) = for all r and p and it follows that / = 0. 



m 3 v„ 
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There remains to consider the term 

C,„, n (k) = f (Z k {\)hxh 2 + Z^hih 2 ) K m ,„ /3 (A, A + h lr A + h 2 ) dA 3 (A, h u h 2 ), 
Jc 3 

where 

Zfc(A) = ^(^G0(A1*). 

«<; 
Observing that supp Z k c supp g and arguing is in the case of A m /H (fc), it is seen that 



f 

Jm, 



Z fc (A) fta ft 2 R m ,nA^ A + h,A + h 2 ) cL4 3 (A, h t/ fc 2 ) -» 0, 



,>3M(S„, 



as ;« — » oo and \n — m%\ < 1. Hence, using [ 12. 5t , we obtain that the asymptotics of C mr „(k) is that 
of C; H/ „(A:) + C^„(fc) where 

C „(*:)= J Z fc (A) fei £ 2 R m ,„, 3 (A, A + hi, A + h 2 ) dA 3 (A, %, fe) = 



= m 3 v 



f Z /C (A) h h 2 (AQ(A) 3 + 0(6,,,)) e'-'AQWC'A-lfcil 2 -^ 2 ) dA3 ( A/ ^ = 

= v m f Z fc (A) ( f (1 + 0(<5„,)) & & e^-l^l 2 -^! 2 dA 2 (4i, &)) dA(A), 

Jc \J|£L<3Mlogm / 

and (likewise) 

(6.6) C,„(/c) = z> m f Z*(A) If (1 + 0(6 m )) |t & e^-^f-l&l 2 ^(^ &) ] dA(A)/ 

Jc \J|4U<3Mlogm / 

where c„, — ♦ 1 as )« — » 00. 

We first claim that C' m n (k) — » when m — > 00 and |n - wit| < 1 for all A: > 2. We will have shown 
that when we can prove that 1/ = where 

L'= f ^V^-^ 2 -^ 2 dA 2 (^ 2 ). 
Jc 2 

To prove this, we pass to polar coordinates and write 

L ' = ^\ J P(r,p)drdp, 
where 

•> 2 7I /^ 2 Tt 



P(r,p)= (rp) 2 e« 6 -V e r P e,(e ^- r2 -P 2 dd d<p. 

Jo Jo 



'0 JO 

Making the change of variables & = 8 — <p and <p = (p, the integral transforms to 

->2n I r>7.n—(p 



Xln I pln-(p \ 

I (rp) 2 e i9 e r P el °dd d<p. 



-<p 
But the inner integral is readily calculated, 

-7.n— <p r -,2n-(p 



Xj.n-(p 
(rp)W*"d$ 



-xrpe'P' 



= 0. 

S=-q> 



This shows that P(r,p) = and consequently U = 0. It follows that C' mn (k) — > as m — > 00 and 
\n -m%\<\ for all k > 2. 

To handle the term C',' n (k), it becomes necessary to calculate 

L"= f U^-\^-^dA 2 {Z, x ,Zi). 
Jc 2 
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Again passing to polar coordinates, we write 

- /-tec /-tea 

L " = ^ I I W(r,p)drdp, 
where 

/-»2ti /-»2ti /-»2ti 

W(r, p) = e~ r2 -P 2 I I (rp) 2 e i(9 -<M e^*" 9 ' d(/> d0 = 2ne~ r2 -P 2 I (rp) 2 e"' 5 e''P e " dS. 
Jo Jo Jo 

We now put z = e lS and use a simple residue argument to get 

m p) = 2Tl(r P )2g " r2 " p2 f ^^ dz = An 2 (rpf e~ r2 -P 2 . 
1 Jt z 

It follows that 
(6.7) 



I (rp) 3 e~ r2 -P 2 dr dp = 1. 
o Jo 



;o ^0 
For ft: = 2 it now follows from J6.7) , | |6.6| | and Lemma [3.4| that 

2 



C" (2) -» - 



f\dg(A) 



dA(A), 



when m — > 00 and |n - ?wt| < 1. On the other hand when /c > 3 we get that 
(6.8) Km C;;„(fc)= fz^A)dA(A) 

m— *c», \n-im\<l J^n 

is pure imaginary, again by Lemma I3.4I In fact this shows that the limit in | |6.8t must vanish, 

because the cumulant C m>ril k(g) is real and all other terms in the expansion (in Lemma |4.3| | but 

C m ,n (k) have already been shown to be real (in fact zero) in the limit when m — > 00 and \n-mx\ < 1 . 

The proofs of all statements are now complete. q.e.d. 

7. Concluding remarks 

We conclude this paper with a series of remarks concerning possible applications and gen- 
eralizations of the main theorem. We also outline an alternative approach to the proof of Th. 



7.1. Non-analytic potentials. Recall that we proved Th. 11.41 assuming that the potential Q is 
real-analytic in some neighbourhood of S % . It is possible to extend this result to more general 
smooth potentials. Assuming that Q is C°°-smooth, one defines the auxiliary functions ip, bo and 
b\ in the expression 

K l m (z,w) = (tnb (z,w) + fci(z,it;))e""' ,(z ' !B) 

as any fixed almost-holomorphic extensions from the anti-diagonal of Q, AQ and ^AlogAQ 
respectively. For example, in the case of i/> this means that ip is well-defined and smooth in 
a neighbourhood of the anti-diagonal in C 2 , and (i) i/> (z, z) = Q(z), (ii) the anti-holomorphic 
derivatives djip vanish to infinite order at each point of the anti-diagonal, i = 1,2, and (iii) 
ip(z, w) — ip(w, z) whenever the expressions make sense. Lenuna [l.2l extends to this more general 
situation; the proof is not very different from the argument in (TJ but it involves some additional 
technical work. The rest of the proof of Th. ll.4l for smooth potentials requires only minor changes. 
As we mentioned earlier, the smoothness (or analyticity) condition is "local" - we need it 
only in some neighborhood of the droplet. In particular, Theorem 11.41 is true for potentials 
Q : C -> R U {+00} of the form 

Q(z) = Q (z) + f log — ^d^(zo), 
Jc \z - Zol 
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where Qo is a smooth function (with sufficient growth at infinity), and y. is a positive, finitely 
supported measure (linear combination of Dirac measures). In this case the droplet S is disjoint 
from supp ,u, and so the "local" smoothness condition holds. (We will need this observation later.) 

7.2. Variational approach. Here we sketch a different, more "physical" proof of our main result, 
Th 11.41 The proof is based on a variational argument well known in the physical literature, see 
e.g. the papers of Wiegmann and Zabrodin. In the rigorous mathematical setting, this method 
was developed by Johansson in the one-dimensional case, see [26|. 

We will use the fact that the estimate (JT7T2J) for K m/n (z, z)e~ m Q^ is uniform when we make 
small smooth perturbations of the potential Q. We will also need some basic facts concerning 
the variation of the droplet under the change of potential (Hele-Shaw theory). Modulo these 
technical issues (see Remark l7.lt the proof of the theorem is rather short. 

To simplify the notation we assume m = n and t = 1 and write K n instead of K n r „, etc. Let 
h : C — » 1R be a bounded smooth function. We denote, for a positive integer n, 

h(z) 
Q„(z) = Q(z) - ±±, 

and we will use "tilde-notation" for various objects defined w.r.t. the weight Q„. Thus K n is the 
kernel function w.r.t. Q„ etc., while the usual notation (K„, etc.) is reserved for the weight Q. 

It is known that, for any K <s S° n X, the coincidence set {Q n = (Qn) }, and therefore the 
perturbed droplet, will contain K in its interior when n is large enough. One can then prove that 

(7.1) K n (z,z)e- n Q»V = nAQ n (z) + ^Alog AQ„(z) + o(l), (n -» oo), 

for z 6 K, and that the o(l)-term is uniform in z. 
Let geCg (Si n X), so we have 

K„(z,z)<r"Q"< 2 > = nAQ(z) - Ah(z) + ^Alog AQ(z) + o(l) 
uniformly for z e supp g. We define 

D',l[g] = E n ({[uct n g). 
If V denotes the Vandermonde determinant, we then have (see (jl .3b and (JTT2J) 
JL fluct„ g-\V\ 2 e - ntmce - Q»dA„ 



D n n [g] 



J c J^| 2 e-" trace " Q "dA H 
f c „ fluct,, g ■ e trace » h | V | 2 e~" trace " QdA„ E„ (fluct„ g ■ e trace » '') 



j^ e trace„ h | V f g -n trace,, Q dAn £„ ( e trace„ fe) 

We now fix a real-valued g and set 

(7.2) h = Atg- JgAQdAY 

where A is a real number, so that 

trace„ h - A fluct,, g. 

We have 

£„ (fluct,, ?-e Afluct »^) 
D «^ = MeAfluct „ g) - F„(A), where F„(A) := logfE^ 3 ^^) . 
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Now from |[7J} we see that 

D h n [g] = [ g(z)K n (z,z)e-"^dA(z) -n f gAQdA = 

= - I A/z-gdA + I gdv + o(l)^> I c?/i • dgdA + J gdv. 

It follows from (|Z2} that 

F' n (A) — > I gdv + — I |Vgj cL4 as n — > oo. 

The last relation can be integrated over A e [0,1]. This is justified by dominated convergence and 
the estimate F" > 0, which is just the Cauchy-Schwarz inequality. It follows that 

log£„ e fluct «* = F„(l) = J F;,(A)dA ^ fgdv + lj\Vg\ 2 dA 

when n — > oo. This means that 

logE„e tfluct "^£e g + £ 2 i> 2 /2 

for all suitable scalars f, which in turn implies Th. 11.41 

Remark 7.1. We have discussed two rather different proofs for our main result Th. 11.41 We remark 
that, in the (most interesting) case when the potential is real-analytic in a neighbourhood of the 
droplet, the theory of asymptotic expansions for the correlation kernel is somewhat simpler and 
cleaner than in the smooth case. In the variational proof we need to make a smooth perturbation 
of the potential, and so we need a discussion of the smooth theory even in cases when the potential 
is real analytic. One would also need to include a further discussion of Hele-Shaw theory to make 
the variational proof complete. We will discuss the variational approach in greater detail in our 
forthcoming paper [2]. 

7.3. Interpretation in terms of Gaussian fields. Denote LI = S{ D X and let < W {U) = W^LI) be 
the completion of CJf (LI) under the Dirichlet inner product 



Jc 



<f,gh= Vf-VgdA 



Let G be the Green's function for LI and denote by £(LI) = W 1,2 (LI) the Hilbert space of distribu- 
tions with inner product 

(pi,p2)&= I I G{z,w)dpi{z)dp 2 {w). 
Ju Ju 

(More accurately, £(LI) is the completion of the space of measures with finite £-norm.) We have 
an isomorphism 

Ay : 1V (U) -> 6(12), 
where Ay = dd is the (Dirichlet) Laplacian. The inverse map is given by the Green potential 
1 



--Au 1 p=U H G where XTJz) = G(z,w)dp(w). 

By a Gaussian field indexed by < Wq{U) we mean an isometry 

r:^Vo(LI)^L 2 (n,P), 

where (Q,P) is some probability space, and T(g) ~ N\0, \\g\\y) for any g 6 < U / o(LI). We now pick 
(Aj)" randomly w.r.t. !!„„ and consider the sequence of random fields (measures) 



r„ = 4 



^ §a ; -iwi-v 

y =1 
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which satisfy 



T n (g)=4\f[uct n g- Jgdvj. 



Thus Th. 11.41 implies that as n — » <x>, the fields r„ converge to a Gaussian field T indexed by 
'H / o(!i)- The precise meaning of the field convergence is convergence of the correlation functions: 

(7.3) e„ (r«(gi) • ■ ■ r„(gfc)) -» <r(gi) • ■ ■ r(#)> 

for all finite collections of test functions {gy} C Cf(U). The right hand side in | |7,3| | is given by the 
Wick's formulas 

<r(g 1 )--r(g V i)> = o 

and 

P 

<r(gi)-r(g 5p )> = £J]< a/ ^> v , 

where the sum is over all partitions of {1, ... , 2p] into p disjoint pairs (4/ /*:)• 
Using the identifications mentioned above, we obtain the following result. 

Proposition 7.2. The random functions 

( " 
h n (z) = 2 ^GfoAyJ-U^z) , 

converge in U to a Gaussian free field with Dirichlet boundary condition, i.e. to a Gaussian field indexed 
byS(U). 

Alternatively, if we pick (Ay) and (A') independently w.r.t. Yl nn then the random functions 

n 

h n (z) = Y J {G(z,A ! )-G(z,A' j )) 

7=1 
converge to a Gaussian free field with Dirichlet boundary condition. 

7.4. Fluctuations near the boundary. In a separate publication [2 J we will prove a version of Th. 
1 1.41 valid for general test functions, which are not necessarily supported in the droplet but just, 
say, of class C^°(C). The proof is based on Ward's identities and Johansson's variational technique 
mentioned above. Here we only settle for stating the result. 

We assume throughout that Q is real-analytic and strictly subharmonic in some neighbourhood 
of the droplet S = Si . One can then prove that the boundary dS is regular, i.e., a finite union of 
real-analytic curves. We will write ds for the arclength measure on dS\ divided by 2n. Denote 

U = S° and U. = C \ S. 

We then have an orthogonal decomposition of the Sobolev space "W = W 1,2 (C), 

f W = < W (U)<B'W(dS)(B'Wo{U>). 

Here ^WoCLT) and ^WoCLT.) are identified with the subspaces of functions which are (quasi- 
everywhere) zero in the complement of LI and IL, respectively, while the subspace "W (dS) consists 
of the functions which are harmonic off dS. The orthogonal projection of *W onto ^{dS), 

f H» f S , 

is just the composition of the restriction operator f (-> /L and the operation of harmonic 
extension to U U U, U {co}. For / 6 "W we also denote by f s the orthogonal projection of / onto 
f W {U)@ f W(dS), 

f s = l s -f + l u ,-f S , 
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in other words, / coincides with f on S and is harmonic and bounded in the complement of 
that set. 

Finally, we write tiuf for the exterior normal derivative of / L and rtu./ the exterior normal 
derivative of / ,, • We can now state the theorem. 

Theorem 7.3. Let f e C^°(C). Then the random variables fluct,,/ on the space (C", Yl n ,n) converge in 
distribution toN[ef,v 2 ), where 



and 

e f 



f/dv + i f n u {f)ds + \ f (/• n u . (log AQ) -n u . (/**)• log AQ)ds. 



Note that the formula for ef becomes very simple in the case of the so called Hele-Shaiv 
potentials , i.e. if AQ = const. > in a neighbourhood of S, then 



(7.4) e f = \ f n u (f)ds 

4 J3S 



In field theoretical terms, Th. [73] means that the random measures 

( " 
^ 6a ; -noi-v 

converge in C to the sum of two independent Gaussian fields - indexed by *Wq ( U) and by *W {dS\ ) 
respectively. While the first one is conformally invariant, the second one is not. 
Alternatively, we can say that the random functions 

p(z;Mi) 



h n (z) = log 



viz; Mi) 

where the p(z; Mi) are the characteristic polynomials of two independent n X n random normal 
matrices My, converge to a free Gaussian field on S with free boundary condition. 

7.5. Large volume limit. Let us take a point zq e vS^nX and assume for simplicity that AQ(zo) = 1. 
Define /.(„ e Prob(C") as the image of Yl„ /tl under the map 

(Ay)^^(V^(Ay-zo))" =1 , 

and think of \x n as a point process in C 

Proposition 7.4. The processes \i n converge to the Ginibre(oo) point process, i.e. to the determinantal 
process with correlation kernel 

K{z,w) = e m -W 1+ W 2)l2 . 

Proof. Assume w.l.o.g. Zq - 0. Then p n are determinantal processes with correlation kernels 

k n (z,w) = -K nn l— ,— ) e -"(QWVi) + QWVS))/2, 
n ' \ V« yn) 
Using the expansion for X„ ; „ in Lemma [i~2l we see that 

k n {z,w) = (AQ(0) + o(l))e n,/ '( z/V "' ir ' /V ")-"( Q(2/V " )+Q(a ' /V " ) ) /2 , 

where the o(l) is uniform for z and w in a fixed compact subset of C Next observe that, up to 
negligible terms, we have 

ty (z, w) = Q(0) + az + aw + bz 2 + bw 2 + zw, 
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for some complex numbers a and b. It follows that 

k n (z w) = (1 + (i)) e i V" Im ( fl ( z - l ")) e iIm ( fc ( z2 ~ n ' 2 )) e zlI '~(l z l 2+ l n; l 2 )/ 2 . 

The first two exponential factors cancel out when we compute the determinants representing 
intensity fc-point functions, which yields the desired result. □ 

7.6. Berezin transform and fluctuations of eigenvalues. We will write 



R*(Ai, . . . , At) = det (K n (A lV A,-)) r» E * Q(A 



for the fc-point intensity function of the ensemble (11.2b with m = n. We will also need the connected 
2-point function 

R 2 f{z,w) = R 2 n (z,w) - Rl(z)Rl(w) = - \K„{z,wf e- n(Q(z)+Q(w)) . 

It is easy to check that 

[ R 2 n ' c (z,zv)dA(w) = -Rliz), 
Jc 

and 

Gov (fluct„/, fluct,, g) = ff(z)g(z)Rl(z)dA(z)+ f f(z)g(w)R 2 n ' c (z,w)dA 2 (z / w). 
Jc Jc 2 

Recall that for a given z, the corresponding Berezin kernel B,, is given by 

Rl(z) R\(z) 

and the Berezin transform is 

Bnf(z)= [ f(w)B<*\w)dA(w). 

Jc 

We may now conclude that 

Cov (fluct,, /, fluct,, g) = f (/(z) - S„/(z)) g(z)Ri(z)dA(z). 
Jc 

On the other hand, Th. I1.4l implies that 

Cov„(fluct„/,fluct„g)-»- I Af(z)g(z)dA(z), (n-»oo), 

Jc 

where /, g 6 C™ (S\ n x). Therefore, 

J(/(z) - S„/(z))Ri,(z)g(z)dA(z) -» - j Af(z)g(z)dA(z). 

Since 

R^nAQ+iAlogAQ(z) + (l) 

on the support of g, we obtain the following asymptotic formula for the Berezin transform. 
Proposition 7.5. Iff 6 C^° (<S° n x), then 

(7-5) S "/ = / + ^ + S 

znsz'de f/ze droplet in the sense of distributions. 
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Berezin's transform has the following probabilistic interpretation. Let us think of the measure 
n„ = n„„ as the law of a point process O,, in C. We will refer to 0„ as the « -point RNM (random 
normal matrix) process associated with potential Q. 

Let us now condition O,, on the event {zo e O,,} and write 0?% for conditional (n - l)-point 
process. Accordingly, we write R k n for the fc-point intensity function of <5„ and R* 1 = RJT the 
fc-point function of ®^\. 

Lemma 7.6. 

(7.6) B< 2 ° > (z)=^(2)-Ri_ 1 (2). 

Proof. Consider small discs D and Do centered at z and Zo with radii e and eo respectively. We 
have 

n„ (jo„ n D * 0j) 



Rj(z ) = lim 



«d-0 £ 2 



and 



It follows that 



n„ ({o„ nD±®}n {<t>„ nD ^ 0}) 

R^(z ,z) = lim 



■ ■-" e 2 ^ 



n„({<t>„nD#0|<t>„nD o *0}) 

R^(z) = lim lim — i! -! !i = 

n„ ({o„ n D * 0} n h„ n D * 0!) 

= lim lim ^ -—!■ [ - - >-?- 

£ ^o fo -o E 2 Un (j$ n n Dq ^ 0jj 

_ R«(Zp,z) _ R i, v p<z >/_ -, 

- Dl , : - K n( z ) - b n \ z >- 



Integrating J7.6I I against test functions we get the following formula, where E n stands for the 
expectation with respect to n„ and E_l with respect to the law of &_l. 

Corollary 7.7. Let z Q eCandf 6 C b (C). Then 

Bnf{zo) = E n (trace,, /) - fig (trace„_i /) . 

The central limit theorem for Berezin transform states that the rescaled (as in the large volume 
limit procedure) Berezin's measures converge to the standard Gaussian distribution in C, see 
|Hl,Th. 2.6. We can now interpret this statement in terms of random eigenvalues. 

Let zo 6 S^P\ X and assume w.l.o.g. that AQ(zo) = 1. Define <J>r°i as a point process in C 
obtained from O^: by dilating all distances to Zo by a factor of yfn as in the previous subsection. 
In other words, we condition O n on the event "zo is an eigenvalue" and rescale the distances. 

Proposition 7.8. The limiting point process of®„, (n — » oo), has the following one-point intensity 
function: 

R 1 ' <2o) (z) = l-e- |z - z ° |2 . 

Proof. Let R '_ ?° ' denote the one-point function of O _l. Similarly, let Rj be the one-point function 
for the process O,,, by which we mean <5„ dilated by a factor of yfn about zq. By Proposition 8.3, 
the point processes $> n converge to Ginibre(oo) ensemble as n —* 00. The one-point function of 
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Ginibre(oo) is R 1 (z) = 1 and its Berezin kernel is B <2o> (z) = e |z 2 °' . Conditioning the equation j7.6\ 
on the event "zq is an eigenvalue", we get 

B^\z) = Rl(z)-E};^(z), 
and sending n — » oo we get the stated formula. n 

7.7. Berezin transform in quasi-classical limit and orthogonal polynomials. As before, let O h 
be the n-point RNM process associated with potential Q. We fix a point zq and condition O n on 
the event {z 6 <5„}. 

Lemma 7.9. T/ze conditional (n - l)-point process <&_* is the RNM process associated with the potential 

Q(z) = Q(z) - -i^ (log |z - z | 2 - Q(z)) . 

Proof. The density of the measure Tl n is given by 

(7.7) p{A lr ...,K) = \ W„(A 1 A„)| 2 e -»(Q(^)+...+Q(A»)) / 

where Z is the normalizing factor (partition function) and V„ the Vandermonde determinant, see 
(11.3b . Setting zo = A„, we have 

p(Aa A„_!, zo) = ?__ |y„_ l(Al A,^)! 2 e^QM^Q^HZP iog|A,- Zo | 

= ^|y„_ l( A 1 A^iv^N^-^M. 

It follows that the density of the conditional point process ^„°\ is 

p(A a A,,-!) = I IV^CAx A„_0l 2 e H^&C*0 + ... + 5.<M, 

z 

where Z is the corresponding normalizing factor. □ 

Let us now assume that the potential Q is real analytic and strictly subharmonic in some 
neighbourhood of the droplet S = Si so that Theorem [73] applies. Denote 

Q„(z) = Q(z) - ^, h(z) := log |z - z | 2 - Q(z). 

i.e. so that Q„ = Q - /z/n. As in Subsection l7.21 for a bounded smooth function / we write 

D„ \f] = E„ (fluct„ /) , D h n [f] = E„ (fluct„ /) , 

where E„ is the expectation with respect to the potential Q n . 

The argument in Section 8.2 shows that the variance part of Theorem [73] is equivalent to the 
statement that 

D n \f\-l%lf\->\<f s ,hh, 

where f s is the orthogonal projection of / onto < Wq( U)®'W(dS) . By Corollary l7. 7l and Lemma [7!9l 
we have 

S„/(zo) = E„ (trace„ /) - E»_i (trace„_i /) 

= ffda + E n (fluct,, /) - E„_! (fluct„_! /) 

= J/da + DJ/J-DjUl/], 
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and therefore 

(7.9) S„/(z ) -» jfdo + (f,h) v , (n -» oo). 

Note that 

(f,h)v = (f,Q S )v-(f,l)v, 
where /(z) = log \z — Zq| and 



(f, Of >v = - JfA&dA = - Jfdo, 



and 

-(f,lh = ff s AldA = f{z ). 



In view of (|7.9t , it follows that 

S„/(zo) -» f(zo). 
Since the function / was arbitrary, we have derived the following result. 

Theorem 7.10. Let zq e C. Then the Berezin measures B„dA converge to the Dime measure at zo if 
z e S\, and to the harmonic measure ofC \ Si evaluated at Zq if Zo £ S\. 

This theorem is also true at Zo = oo, in which case it has the following form. 

Theorem 7.11. Let P„ be the n-th orthonormal polynomial with respect to the measure e~ n ®dA in C. 
Then the probability measures 

\P n \ 2 e~ nQ dA 
converge to the harmonic measure ofC \ Si evaluated at oo. 

Proof. We need to compute the limit of the Berezin kernel B„ (z) as zq — » oo. By Lemma [7. 61 we 
have 

B^(z)=R 1 n (z)-Rl 1 (z), 

where R}, and R 1 , are the 1-point functions of O,, and <S>_ 1 respectively. Since O,, is the n-point 
RNM process associated with potential Q, we have 

Rl^Ptfe-^. 

k=0 

On the other hand, by Lemma |7!9l O _. is the (n - l)-point RNM process associated with the 
potential 

2 



Q< Z o> (2) = » Q( Z ) + J- log f _W 

«-l^ w n-1 & \|z-z | 2 



(Here we added a constant term to the potential Q in Lemma 17.91 this clearly didn't affect the 
point process.) Since 

Q <2o> (z) -» Q(z) := -ArQ(z) as z -» oo, 
n — 1 

we have 

limR^^yiP,! 2 ^"-^, 

k=0 

where {P/ c } are orthonormal polynomials with respect to the weight 

e -(«-l)Q _ e -nQ_ 
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Since the weight is the same for the polynomials {P;t} and {pt} we have 
B^ = £ \P k \ 2 e-' lQ ~ £ \P k \ 2 e-' lQ = |P„-il 2 e~"Q. 

k=0 k=0 

Combining this with Theorem 17,101 we conclude the proof. □ 

7.8. Further remarks on the cumulant method. We here continue our discussion of the cumulant 
method (Sect[TJ and compare our result with some other related work using this method. 

In |35J, Soshnikov studied linear statistics of the form trace,, g n - E(trace„ gn) where g„(t) = 
g(L n t) and L n is a fixed sequence with L„ — > co,L n /n — » 0. The expectation is here understood with 
respect to the classical Weyl measure on [—n, n) n , i.e., we are considering the Gaussian unitary 
ensemble; g : M — » 1R is a test function in the Schwarz class. 

In 1 35 J, asymptotic normality is proved for these linear statistics using the cumulant method 
applied to the sine-kernel, i.e. the explicit correlation kernel in that case. The asymptotic 
variance of trace„ g„ turns out to be finite and independent of the particular sequence L„; it equals 

The method in ll35l does however not allow to draw conclusions about the case L„ « 1; the 
assumption L„ — » oo is used in the proof of Theorem 1 (p. 1357), where limits of certain Riemann 
sums are identified. 

We also want to mention the short proof of asymptotic normality due to Costin and Lebowitz 
lll6l . In the situation of [16J, one considers certain linear statistics which have infinite asymptotic 
variance. This infiniteness of the variance is then used to show decay of the cumulants of the 
corresponding normalized variables. (Thus the method in |16| necessarily breaks down in our 
situation, when the variance tends to a finite limit.) 

The cumulant method has also been used in the theory of Gaussian analytic functions, see 1311 . 
In this case, asymptotic normality was obtained for linear statistics whose variances converge to 
zero. In [33], the result was generalized to a setting of zeros of random holomorphic sections of 
high powers of a positive Hermitian line bundle over a Kahler manifold. (Cf. the book l23l for 
further developments in the theory of Gaussian analytic functions.) 

Acknowledgements. We are grateful to Alexei Borodin, Kurt Johansson and Paul Wiegmann for 
help and useful discussions. 
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